NTUTH-03-505C 
May 2003 



Light quark masses, chiral condensate and 
quark-gluon condensate in quenched lattice 
QCD with exact chiral symmetry 

Ting-Wai Chiu and Tung-Han Hsieh 

Department of Physics, National Taiwan University 
Taipei, Taiwan 106, Taiwan. 
E-mail : twchiu@phys.ntu.edu.tw 



Abstract 

We determine several quantities in quenched lattice QCD with exact 
chiral symmetry. For 100 gauge configurations generated with Wilson 
gauge action at f3 = 6.0 on the 16 3 x 32 lattice, we compute quenched 
quark propagators for 13 bare quark masses. The pion decay constant 
is extracted from the pion propagator, and from which the inverse lat- 
tice spacing is determined to be a -1 = 1.979(6) GeV. The parameters 
(8, A, B) in the pseudoscalar meson mass formula l)14|) in quenched chiral 
perturbation theory (q;\PT) to one-loop order are determined. Further, 
we measure the index (topological) susceptibility of these 100 gauge 
configurations, xt = (175 ±6 MeV) 4 , from which we obtain an estimate 
of the mass of rf in q^PT, and the coefficient of quenched chiral loga- 
rithm, both in good agreement with the values determined from the pion 
masses, as well as with the theoretical estimates. With our values of 
5, A, B, the experimental inputs of pion and kaon masses, and the pion 
decay constant, we determine the light quark masses: m u ^ = 4.1 ± 0.3 
MeV, and m s = 92 ± 9 MeV, in the MS scheme at scale fi = 2 GeV. 
Also, we determine the quark condensate (qq) = —(250 ± 3 MeV) 3 , and 
the quark-gluon condensate giq&^F^q) = — (434±4 MeV) 5 , in the MS 
scheme at scale 2 GeV. 
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1 Introduction 



Exact chiral symmetry on the lattice was pioneered by Kaplan il] with his pro- 
posal of domain-wall fermions (DWF) on the 5-dimensional lattice, in which 
the fifth dimension (internal flavor space) is discretized with N s sites (flavors) 
and lattice spacing a^. Although the initial motivation was to provide a non- 
perturbative formulation of chiral gauge theories, the idea turns out to be 
natural and well-defined for vector-like gauge theories (e.g. QCD), with quark 
fields constructed from the boundary modes with open boundary conditions 
Soon after Kaplan proposed DWF for chiral gauge theories, Narayanan 
and Neuberger j^j observed that the chiral determinant can be written as the 
inner-product ("overlap") of two fermionic many body states of two bilinear 
Hamiltonians. For vector gauge theories like QCD, the fermion determinant 
is the product of a complex conjugate pair of chiral determinants, thus it is 
gauge invariant, real, non- negative, and the corresponding lattice Dirac opera- 
tor for massless quarks (i.e., the overlap Dirac operator [4 ) can be represented 
by a finite matrix of fixed shape regardless of the topology of the background 
gauge field, without undesired doubling or any fine-tuning. Mathematically, 
the overlap Dirac operator is exactly equal to the effective 4D lattice Dirac 
operator (for internal fermions dressed with pseudofermions) of DWF (with 
m q = 0) in the limit N s — > oo followed by as — > 0, 

D = m ( 1 +75-^= ) , H w = j 5 D w (1) 



where D w is the standard Wilson Dirac operator plus a negative parameter 
— m (0 < m < 2). 

However, it has been shown jH] that the conventional DWF with any finite 
N s does not preserve the chiral symmetry optimally for any gauge background. 

Recently, one of us (T.W.C.) has proposed a new DWF action E] such 
that the quark propagator preserves the chiral symmetry optimally for any N s 
and background gauge field. Further, its effective 4D lattice Dirac operator 
(with any N s and m q ) for internal fermion loops is exponentially-local for 
sufficiently smooth gauge backgrounds jjj. In the optimal DWF, the quark 
fields are constructed from the boundary modes at s = and s = N s + 1 6J, 

q( x ) = yfi [P-i/>{x, 0) + P + ^(x, N s + 1)] , r = — - (2) 

ZTHq 

q(x) = [$(x, 0)P+ + $(x, N s + 1)P_] , (3) 

and the generating functional for n-point Green's function of the quark fields 
has been derived [Hj, 

_ J[dU]e- A «det[(D c + m q )(l + rD c )- 1 } exp { J(D C + m q )- l j} 
Z ^ JJ ~ f[dU]e- A °det[(D c + m q )(l + rDc)' 1 ] 



where A g is the action of the gauge fields, m q is the bare quark mass, J and 
J are the Grassman sources of q and q respectively, and 



rD c 
S pt 



1 + l5S0.pt 

1 - 75S opt 

1 H w R%- l > n) {Hl), 



N s = 2n + 1, 
7V e = 2n. 



(5) 
(6) 



Here Rz{H^) is the Zolotarev optimal rational polynomial [Hj for the inverse 
square root of H 2 , 
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where the coefficients d , d > Q and are expressed in terms of elliptic functions 
[S] with arguments depending only on N s and X max /X min (X max and A mi „ are 
the maximum and the minimum of the the eigenvalues of |i?u,|). 

From Pjl. the effective 4D lattice Dirac operator for the fermion determi- 
nant is 

D(m q ) = (D c + m 9 )(l + rD^ 1 = m q + (m - m q /2) [l + l5 H w R z {Hl)] , (9) 
and the quark propagator in background gauge field is 

5 2 Z[J, J] 



(q(x)q(y)} 



j=j=o 



5J(x)5J(y) 

(D c + m,)-; = (1 - rm q )- l [D- x l{m q ) - rS x J (10) 



In practice, we have two ways to evaluate the quark propagator (fTU|) in 
background gauge field: 

(i) To solve the linear system of the 5D optimal DWF operator, as outlined in 
Ref. 0; 

(ii) To solve D~^ y {m q ) from the system 



D{m q )Y 



m, 



+ (m - m q /2) (l + l5 H w R z (H 2 u 



Y = 1 



'IV. 



with nested conjugate gradient |9:, and then substitute the solution vector Y 
into (HOD- 
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Since either (i) or (ii) yields exactly the same quark propagator, in prin- 
ciple, it does not matter which linear system one actually solves. However, 
in practice, one should pick the most efficient scheme for one's computational 
system (hardware and software). 

For our present software and hardware (a Linux PC cluster of 42 nodes [TU]). 
it has been geared to the scheme (ii), and it attains the maximum efficiency 
with Neuberger's 2-pass algorithm |TJ for the inner conjugate gradient loop of 
jUJ). So, in this paper, we use the scheme (ii) to compute the quark propagator, 
with the quark fields 0-© defined by the boundary modes in the optimal 
DWF. A study of Neuberger's 2-pass algorithm is presented in Ref. [T9"] . 

For our computations in this paper, we fix mo = 1.3, and project out 20 low- 
lying eigenmodes and two high- lying modes of \H W \. With N s = 32, i.e., 16-th 
order Zolotarev rational polynomial R^ 5 ' 16 \h^) in (|TT|) . we impose stopping 
criteria 2 x 10~ 12 and 10 -11 for the inner and outer conjugate gradient loops 
respectively. Then the chiral symmetry breaking due to finite N s is always less 
than 10" 13 , 



a 



opt i 



< 1(T 1J , (12) 



for every iteration of the nested conjugate gradient. 

The outline of this paper is as follows. In Section 2, we determine the 
parameters (5, A, B) in the pseudoscalar meson mass formula (fT4^) in quenched 
chiral perturbation theory (q;\PT) to one-loop order, by computing quenched 
quark propagators for 13 bare quark masses, and for 100 gauge configurations 
generated with Wilson gauge action at (3 = 6.0 on the 16 3 x 32 lattice. The 
pion decay constant is extracted from the pion propagator, and from which 
the lattice spacing is determined to be a = 0.0997(3) fm = 0.5053(15) GeV" 1 . 
Also, we measure the index (topological) susceptibility of these 100 gauge 
configurations, Xt — (175 ± 6 MeV) 4 , which provides an estimate of m v / = 
813 ± 51 MeV in q;\PT, and the coefficient of quenched chiral logarithm 5 = 
0.16(2), which is consistent with 5 = 0.164(13) determined from the pion 
masses. In Section 3, with our determined values of 5, A, B, the experimental 
inputs of pion and kaon masses, and the pion decay constant, we determine the 
light quark masses: m U ( i = 4.1 ± 0.3 MeV, and m s = 92 ± 9 MeV, in the MS 
scheme at scale fi = 2 GeV. In Section 4, we determine the quark condensate 
(qq) = —(250 ± 3 MeV) 3 , and the quark-gluon condensate giqcr^F^q) = 
-(434 ± 4 MeV) 5 , in the MS scheme at scale 2 GeV. 
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2 Determination of the Parameters in Quenched 
Chiral Perturbation Theory 

In quenched chiral perturbation theory (q%PT) [T3*l IT4*j the pion mass to one- 
loop order reads 

ml = Cm q {l - 5[\n(Cm q /A 2 x ) + 1]} + Bm\ (13) 

where m q denotes the bare (u and d) quark mass, A x is the chiral cutoff, C 
and B are parameters, and 5 is the coefficient of the quenched chiral loga- 
rithm. Further, summing all diagrams in which each internal rj propagator 
is decorated with any number of "bubble sums" (i.e., summing the leading 
[5\n(Cm q /A^)} n terms to all orders), one obtains ^3] 

ml = Amf~ s + Bm\ (14) 

where 



A 



A 2 x 5 



e 



1+5 

(15) 



Obviously, the value of 5 is independent of the chiral cutoff A x (which is usually 
taken to be I^Txf^ with f w ~ 132 MeV). 

In this section, we determine 5 by fitting our data of ml to (j!3j) and (fT4*j) 
respectively, and find that these two values of 5 are in good agreement with 
each other. 

Theoretically, 5 can be estimated as [TB*] 



2 

m 

S =^wr (16) 

where /„- denotes the pion decay constant, Nf denotes the number of light 
quark flavors, and my denotes the rj mass in qxPT, which can be estimated 

as 



m 2 , + m\ - 1m\ = 853 MeV , (17) 

with experimental values of meson masses 1 : m q i = 958 MeV, m v = 547 MeV, 
and m K = 495 MeV. For U = 132 MeV, N f = 3, and my = 853 MeV, (flfi) 
gives 

5 ~ 0.176. (18) 



^^Here we distinguish the physical rj' mass (my), from the 77' mass (my) in q%PT. 
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Evidently, if one can extract 5 from the data of rri^, then m v > in qxPT can 
be determined by (fTBj) . 

Besides from the data of m 2 , one can also obtain S via by extracting 
rrirf from the propagator of the disconnected hairpin diagram. However, to 
compute the propagator of the hairpin is very tedious. 

Fortunately, with the realization of exact chiral symmetry on the lattice, 
the quark propagator coupling to 7/ is (-D c + m g ) _1 , thus only the zero modes of 
D = Dcll+rDc)' 1 can contribute to the hairpin diagram, regardless of the bare 
quark mass m q . Therefore one can derive an exact relation ^lEUE! between 
the rj' mass in qxPT and the index susceptibility of D, without computing the 
hairpin diagram at all. Explicitly, this exact relation reads as 

, 9 AN f ((n + -n) 2 ) AN f , , 

where N is the total number of sites, and Xt = (( n +~ n -) 2 ) /N is the index sus- 
ceptibility of D in the quenched approximation. Then (|T^|) and (|19|) together 
gives 

6 = N (20) 

A salient feature of ()2U|) is that 5 can be determined at finite lattice spacing 
a, by measuring the index (susceptibility) of D, and with f n a extracted from 
the pion propagator. 

Now it is clear that, in order to confirm the presence of quenched chiral 
logarithm in lattice QCD, one needs to check whether the coefficient S ob- 
tained by fitting (fTlj) to the data of m 2 , agrees with that (J2Uj) from the index 
susceptibility. This is a requirement for the consistency of the theory, since 
the quenched chiral logarithm in (fT3j) is due to the 77' loop coupling to the 
pion propagator through the mass term (in the chiral lagrangian), thus 6 {m^) 
must be the same in both cases. We regard this consistency requirement as 
a basic criterion for lattice QCD (with any fermion scheme) to realize QCD 
chiral dynamics in continuum. 

Recently, we have determined 5 = 0.203(14), from the data of m 2 , as 
well as 5 = 0.197(27) from the index susceptibility, for the 8 3 x 24 lattice at 
(3 = 5.8 [IH]. Their good agreement suggests that lattice QCD with exact 
chiral symmetry indeed realizes QCD chiral dynamics in the continuum. 

In this paper, we extend our studies to a lattice of larger volume and 
smaller lattice spacing, 16 3 x 32 at (3 = 6.0. Details of our computational 
scheme have been given in Ref. [T^], except the inner conjugate gradient loop 
is now iterated with Neuberger's 2-pass algorithm [11 J which not only provides 
very high precision of chiral symmetry with a fixed amount of memory, but 
also increases (~ 20%) the speed of our computational system (a Linux PC 
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cluster of 42 nodes at the speed ~ 50 Gflops [TUJ ) . A study of Neuberger's 
2-pass algorithm is presented in Ref. [T2"] . 

We retrieved 100 gauge configurations from the repository at the Gauge 
Connection (http://qcd.nersc.gov), which were generated with Wilson gauge 
action at f3 = 6.0 on the 16 a x 32 lattice [TH] - For each configuration, quark 
propagators are computed for 13 bare quark masses. 

Then the pion propagator 

M 5 (x, t; 0, 0) = tr{ 75 (L> c + m,)" 1 ^, x) l5 {D c + m q )'\x, 0)} , (21) 

(where the trace runs over the Dirac and color space), and its time correlation 
function 

G(t)=^M 5 (x,t;0,0) , (22) 
are obtained, and (G(t)) is fitted by the usual formula 

2m 7r a 

to extract the pion mass m^a and the pion decay constant 

f n a = 2m q a 



In Figs. Qand El we plot the pion mass square (m^a) 2 and decay constant 
f n a versus bare quark mass m q a, respectively. 
The data of f w a (see Fig. |2J) can be fitted by 

Ua = 0.0667(2) + 0.2218(20) x (m q a) 

with x 2 /d.o.f.=0.95. Thus taking f n a at m q a = equal to 0.132 GeV times 
the lattice spacing a, we can determine the lattice spacing a and its inverse, 

1 32 

a' 1 = — — GeV = 1.979(6) GeV , 
Jo 

a = 0.0997(3) fm . (23) 

Thus the size of the lattice is ~ (1.6 fm) 3 x 3.2 fm ~ (1.9 fm) 4 . Since the 
smallest pion mass is 412 MeV, the lattice size is ~ (3.3) 3 x 6.6, in units 
of the Compton wavelength (~ 0.48 fm) of the smallest pion mass. For the 
largest pion mass (1052 MeV), the lattice spacing is about 0.5 of its Compton 
wavelength. 

Next, rewriting ()14|) as 

mla 2 = A 1 (m q a)~tTs + B(m q a) 2 , A 1 = AgTW 
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0.35 




0.00 0.05 0.10 0.15 0.20 0.25 

m q a 

Figure 1: The pion mass square (m^a) 2 versus the bare quark mass m q a. The 
solid line is the fit of Eq. (gjj) . 



and fitting to the data of (see Fig. |TJ), we obtain 

S = 0.164(13) (24) 
A 1 = 0.827(25) (25) 
B = 1.886(87) (26) 

with x 2 /d.o.f.=0.54. 

Even though the quenched chiral logarithm may not be easily detected in 
the graph of (m n a) 2 versus m q a, it can be unveiled by plotting (m n a) 2 / (m q a) 
versus m q a, as shown in Fig. El Further, plotting log[(m 7r a) 2 /(m g a) — B{m q a)] 
versus log(m g a), the presence of quenched chiral logarithm is evident, as shown 
in Fig. |U 

Evidently, the coefficient of quenched chiral logarithm 5 = 0.164(13) is in 
good agreement with the theoretical estimate 5 = 0.176 [Eq. (JTHJ)] in q%PT. 
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Figure 2: The pion decay constant f n a versus the bare quark mass m q a. The 
solid line is the linear fit. 



Also, from ([16)1 . we obtain an estimate of the mass in qxPT, 

Try = 823 ± 33 MeV > ( 27 ) 

which is in good agreement with the theoretical estimate ([17)1 . 

At this point, it may be instructive to find out what is the value of 5 if 
one fits (fT3|) to the data of (m^a) 2 , with a fixed A x , say A x a = 2v / 27r/ 7r a = 
2y/2n x 0.0667. The results (with x 2 /d.o.f.=0.87) are 

5 = 0.168(17), (28) 
Ca = 1.093(20), (29) 
B = 2.092(125), (30) 

where the value of 5 (|28|) is in good agreement with (j2H), and other parameters 
are also consistent with those determined with (jl4j) . 

Note that, in Fig. El the minimum of (m^a) 2 / (m q a) occurs at m q a ~ 0.09, 
which corresponds to m q ~ 180 MeV for a ~ 0.0997 fm. This is consistent 
with the data of 10 3 x 24, and 12 3 x 24 lattices at 3 = 5.8 [T3j. where the 
minima of (m n a) 2 / (m q a) occur at m q ~ 190 and 217 MeV respectively. That 
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Figure 3: (m^a) 2 / (m q a) versus the bare quark mass m q a. 



is, the minima of (m T a) 2 / \m q a) for these three lattices occur in the range of 
m q ~ 200 ± 20 MeV. This suggests that one should be able to observe the 
quenched chiral logarithm behavior (in the plot of (m^a) 2 / (m q a) versus m q a) 
for m q < 200 MeV (i.e., m w < 600 MeV). 

However, in a recent study of chiral logs in quenched lattice QCD with 
Iwasaki gauge action [20^, the location of the minimum of (m^a) 2 / (m q a) occurs 
at m q a ~ 0.13 with a ~ 0.2 fm, which corresponds to m q ~ 130 MeV (i.e., 

~ 480 MeV). This seems to suggest that the location of the minimum of 
(m^a) 2 / (m q a) depends on the gauge action as well as the lattice spacing a. 
For a recent review of quenched chiral logarithms of various quark actions and 
gauge actions, see Refs. [2U 1221 and the references therein. 

In order to provide a self-consistency check of the coefficient of quenched 
chiral logarithm (5) determined from the data of m 2 , we measure the index 
susceptibility of D by the spectral flow method, and then determine 5 via (J20j) . 

The distribution of the indices of D for these 100 gauge configurations is 
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legume) 2 / (m q a) - B (m q a)] 
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0.1 
log(m q a) 



Figure 4: The extraction of the quenched chiral logarithm by plotting 
\og[(m n a) 2 1 (m q a) — B{m q a)] versus log(m g a). The slope of the fitted straight 
line is equal to —5/ (1 + 5). 



listed in Tabled So the index (topological) susceptibility is 

((?!+ -n„) 2 ) 



4 

a Xt 



Xt 



N 

175 ± 6 MeV) 4 



6.03(75) x 10" 



(31) 
(32) 



where N is the total number of sites. 

Now substituting the index susceptibility JSJ), the lattice spacing 
f n = 132 MeV, and Nf = 3, into the exact relation ()19j) . we obtain the rj mass 



Try = 813 ±51 MeV , 



(33) 



which agrees with the theoretical estimate (fT7jl . Next we substitute (|3*T|) and 
/ T a = 0.0667(2) into (HI), and get 



5 = 0.16 ±0.02 



(34) 
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n + — n_ 


number of configurations 


5 


6 


4 


5 


3 


6 


2 


8 


1 


14 





10 


-1 


15 


-2 


15 


-3 


10 


-4 


4 


-5 


3 


-6 


3 


-7 


1 



Table 1: The distribution of the indices of D for 100 gauge configurations at 
(3 = 6.0 on the 16 3 x 32 lattice. 

which is in good agreement with the value (|28|) determined from the data of 
m 2 , as well as with the theoretical estimate (fTHj) . 

Note that a finite lattice must impose a lower bound for the pion mass, as 
well as the corresponding one for the bare quark mass. If one decreases the 
bare quark mass beyond its lower bound, then the resulting pion mass would 
be flattening (or increasing) rather than decreasing. This is essentially due to 
the finite volume effects of the zero modes of the lattice Dirac operator, which 
is proportional to \n + — n_|/ \m 2 q a 2 N) as m q a — > 0. Thus, only in the infinite 
volume limit (N — > oo), one can obtain zero pion mass with zero bare quark 
mass. Nevertheless, for a finite lattice, the finite volume effects of the zero 
modes can be suppressed if the bare quark mass m q a 3> \J\n + — n_\/N . For 
a lattice of size 16 3 x 32 at (3 = 6.0, we find that m q a > 0.03 is sufficient to 
suppress the finite volume effects of the zero modes. This can be verified by 
extracting the "pion" mass m^a from the correlation function of the "pion" 
propagator with 7475 coupling 

tr[747 5 (D c + m 9 )~ 1 (0, 2)7475(1^ + m q )~ 1 (x, 0)] , 

which is void of the contributions of topological zero modes, and then compare 
m^a with m n a. In Fig. m^a is plotted against m^a, for the same m q a. The 
good agreement between 1714a and m n a for each of the 13 bare quark masses 
provides an estimate that the error due to the finite volume effect of the zero 
modes is less than 5%, thus is well under control. 

In passing, we present our data of the vector meson mass. In Fig. El we 
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0.6 
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0.4 - 



0.3 



0.2 - 




Figure 5: The "pion" mass m 4 a versus the pion mass m n a. The solid line is 
m 4 = m^. 



plot m p a versus m^a. At the smallest bare quark mass m q a = 0.03 in our 
data, the ratio m^jm p is 0.435(8). Fitting our data to the vector meson mass 
formula in qxPT, 



m p a = C + 5C 1 / 2 (m 7T a) + Ci(m 7r a) 2 + C 3 / 2 (m 7r a) 3 , (35) 

we obtain 

C = 0.474(6) 

SCi/2 = -0.223(56) 

Ci = 1.370(158) 

C 3/2 = -0.542(141) 

with xVd.o.f. = 0.92. 

On the other hand, if we plot m p a versus (m n a) 2 , as shown in Fig. [?l then 
the data seems to be well fitted by 

m p a = 0.495(4) + 4.57(23) x (m n a) 2 (36) 

with xVd.o.f. = 0.98. 
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With our present data, it seems to be unlikely to rule out either one of 
above two possibilities. We hope to return to this problem elsewhere. 

Note that even though one has no technical difficulties to use smaller quark 
masses to push m^/rrip down to smaller values, say, 0.2, one suspects that the 
finite size effects (especially for m n ) might be too large to be meaningful. 
Further, for consistency, all quantities in this paper are determined with the 
same set of quark propagators computed for 13 bare quark masses ranging 
from m q a = 0.03 to 0.2. 

0.75 j 
0.70 - 
0.65 - 

m p a 

0.60 - 
0.55 - 
0.50 - 
0.45 -■ 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

m n a 

Figure 6: The vector meson mass m p a versus the pion mass m n a for 13 bare 
quark masses. The solid line is the fit to Eq. (J35|) . 



3 Light Quark Masses 

In the standard model, the quark masses are fundamental parameters which 
have to be determined from high energy experiments. However, they cannot 
be measured directly since quarks are confined inside hardrons, unlike an iso- 
lated electron whose mass and charge both can be measured directly from its 
responses in electric and magnetic fields. Therefore, the quark masses can only 
be determined by comparing a theoretical calculation of physical observables 
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0.60 



0.55 - 



0.50 



0.45 



0.40 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Figure 7: The vector meson mass m p a versus the pion mass square (m^a) 2 for 
13 bare quark masses. The solid line is the linear fit, Eq. (|3l)J). 



(e.g., hadron masses) with the experimental values. Evidently, for any field 
theoretic calculation, the quark masses depend on the regularization, as well 
as the renormalization scheme and scale. 

One of the objectives of lattice QCD is to compute the hadron masses non- 
perturbatively from the first principle, and from which the quark masses are 
determined. However, the performance of the present generation of comput- 
ers is still quite remote from what is required for computing the light hadron 
masses at the physical (e.g., m w ~ 140 MeV) scale, on a lattice with enough 
sites in each direction such that the discretization errors as well as the finite 
volume effects are both negligible comparing to the statistical ones. Never- 
theless, even with lattices of moderate sizes, lattice QCD can determine the 
values of the parameters in the hadron mass formulas of the (quenched) chiral 
perturbation theory. Then one can use these formulas to evaluate the hadron 
masses at the physical scale, as well as to determine the quark masses. 

In quenched chiral perturbation theory JH1E|, the ratio of pion and kaon 
masses to one-loop order reads 



m 



K 



mi 



m + m s 
2m 



1 + 5 



1 - 



m 



m. 



-In (^i 

m \ m 



(37) 
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where m denotes u and d quark bare masses in the isospin limit (m u = ma = 
m), m s the s quark bare mass, 5 the coefficient of the quenched chiral loga- 
rithm. 

With experimental values of pion and kaon masses as inputs to (J38j) and 
()37j) . one can determine the ratio of light quark bare masses m s /m, but the 
absolute scale cannot be fixed by quenched chiral perturbation theory. At this 
point, lattice QCD plays the important role to fix the values of the parameters 
A, B and 5 in the pion mass formula (fT4*J) . 

ml = Am^Ts + Bm 2 , (38) 

by measuring the pion mass versus the bare quark mass. Then the light quark 
masses m and m s can be determined with experimental inputs of m n , mx, and 
f n (to fix the lattice spacing a). 

In Section 2, we have determined the lattice spacing a, and the parameters 
5, Ai and B (J23"|) -(f2l) |) . by computing the quenched quark propagators for 100 
gauge configurations generated with the Wilson gauge action at ft — 6.0 on 
the 16 3 x 32 lattice. Now inserting 5 (|2*Hjl and experimental values of meson 
masses {jxik = 495 MeV and m n = 135 MeV) into Eq. (jSTj), we obtain the 
quark mass ratio 

Tfl 

— = 22.58(23) , (39) 
m 

comparing with the ratio at the zeroth order (5 = 0) 

fUM = 25.89 (40) 
\m J o 

From l|BBj). with values of a, 5, A x and B in ((23)-®, and m n = 135 MeV, 
we obtain 



m = 4.7 ±0.5 MeV 
which is inserted into (J3T?|) to yield 



m s 



107 ± 11 MeV 



(41) 



(42) 



In order to transcribe above results (|4*T|) - ((4*2*|) to their corresponding values 
in the usual renormalization scheme MS in high energy phenomenology, one 
needs to compute the lattice renormalization constant Z m = Zj 1 , where Z s 
is the renormalization constant for ipip. In general, Z m should be determined 
non-perturbatively. In this paper, we use the one loop perturbation formula 



1 + 



47T 2 



ln(aV) + 6 - 7722 



(mo = 1.30) 



(43) 
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to obtain an estimate of Z s . At f3 = 6.0, a = 0.0997(3) fm, and /i = 2 GeV, 
()43j) gives Z s = 1.172(1). Thus the quark masses in (|41j ) -([42 p are transcribed 
to 



m*f(2 GeV) = 4.1 ± 0.3 ± 1.0 MeV (44) 
m^ S (2 GeV) = 92 ± 9 ± 16 MeV (45) 

where the first errors are statistical errors, and the second errors are our es- 
timate of systematic errors due to finite volume and discretization effects, by 
comparing our present results to our earlier ones [22] with smaller volume 
and larger lattice spacing. Evidently, the light quark masses in (jUJ-tJUJ) are 
in good agreement with the current lattice world average [22], as well as the 
strange quark mass determined in Ref. |2"E] . 



4 Chiral Condensate and Quark-Gluon Con- 
densate 

4.1 Chiral Condensate 

The chiral condensate in lattice QCD with the optimal DWF can be expressed 
as 

1 S S 

— (qq) = lim lim — — - - , Z[J : J] 
™ -^o^oo n ^ SJ(x) 5J(x) 



iq- 



J=J=0 



1 J[dU] e-** detD(m q ) J2 X tr(£> c + m q )~} x 

= lim lim — 

m^o^oon J[dU] e~ A 3 det D(m q ) 

= E(JV» (46) 

where Q = Na 4 is the volume of the 4D lattice, Nf is the number of quark 
flavors, Z[J,J] is the generating functional in (j3J), and the Dirac, color, and 
flavor indices have been suppressed. 

The subtlety of ()46|) lies in its two limiting processes. If the order of these 
two limiting processes is exchanged, then (qq) must be zero since sponta- 
neous chiral symmetry could not occur on a finite lattice for exactly massless 
fermions. Thus it seems to be unlikely to be able to extract the chiral conden- 
sate of QCD on a finite lattice, not to mention in the quenched approxima- 
tion 2 . Nevertheless, using chiral perturbation theory or random matrix theory, 
the chiral condensate at finite volume and quark mass can be computed in the 
continuum, both for unquenched and quenched QCD. In the quenched approx- 
imation, the chiral condensate with fixed m q , Q, and Q (topological charge) 

2 Recall that the 77' loop in the quark propagator leads to the quenched chiral logarithm 
in the chiral condensate. 
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reads |2Zj 

- (qq)m,n,Q 



^IQlO^IQlO) + I \Q\+i(z)K\ Qhl (z) 



(47) 



where £ = E(iVf = 0), z — m g f2£, and Jq and Kg are modified Bessel 
functions. Thus, the chiral condensate can be extracted by measuring (qq) in 
different topological sectors at different m q and Q. Alternatively, a plausible 
way to extract the chiral condensate on a finite lattice (with large volume Q) 
is to measure tr(D c + m q )~\, i.e., 



- (qq) 



f[dU]e- A ° ti(D c + m q )-l 
J[dU] e- A " 



(48) 



for a set of small quark masses m q satisfying z ^> 1. Then the chiral condensate 
at m q = can be obtained by linear extrapolation, since tr(D c + m g )~J, behaves 
like a linear function of m q for small m„. This ansatz can be verified in practice. 



a <q q> 




o.oo 



0.05 



0.10 



0.15 



0.20 



0.25 



m q a 



Figure 8: The chiral condensate a 3 (qq) versus the bare quark mass m q a. The 
solid line is the linear fit (|49jl . 



In the following, we measure (qq) in the quenched approximation, with 
one fixed site at the origin (0,0), for 100 gauge configurations generated with 
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Wilson gauge action at f3 = 6.0 on the 16 3 x 32 lattice. The results are plotted 
in Fig. IHlfor 13 bare quark masses, and they are very well fitted by the straight 
line 

- a 3 (qq) = 1.73(3) x 1(T 3 + 0.242(0) x (m q a) . (49) 

Now we take the value 1.73(3) x 10~ 3 at m q a = to be the chiral condensate, 
then we obtain 

(qq) = -0.0134(2) GeV 3 (50) 

with the inverse lattice spacing a -1 (J23j) determined in Section 2. Using the 
one loop renormalization Z s (|4*3*|) . we transcribe the chiral condensate 1)50)1 to 
MS at scale fx = 2 GeV, 

(qq) MS (2 GeV) = -(250 ± 3 MeV) 3 , (51) 

which is in good agreement with the phenomenological estimate, as well as 
other recent lattice results 123 • Nevertheless, a complete picture of the 
underlying mechanism of spontaneous symmetry breaking in QCD can only 
be unveiled in the context of dynamical quarks. 



4.2 Quark-Gluon Condensate 

The quark-gluon condensate giqcr^F^q) is one of the nonperturbative features 
of the QCD vacuum. It is one of the parameters in the framework of QCD 
sum rule to determine the low-energy properties of QCD |3T?1 13T]. 

The first measurement of the quark-gluon condensate in lattice QCD 
was performed more than 15 years ago, with the staggered fermion, for only 
5 gauge configurations. A recent measurement jHS| was also carried out with 
the staggered fermion, at each site along the body diagonal of the 16 4 lattice, 
for 100 gauge configurations generated with Wilson gauge action at f3 = 6.0. 
In Ref. jSSI, with 1600 measurements for each quark mass (of total 3 quark 
masses), the quark condensate and quark-qluon condensate at m q = were 
obtained by linear extrapolation, however, the ratio g(qc^ v F^ u q) / ' (qq) turns 
out to be much larger than the value in the QCD sum rule as well as other 
phenomenological estimates. 

In general, the quark-gluon condensate in lattice QCD with the optimal 
DWF can be expressed as 

-g^nvF^q) 

1 5 5 
= lim lim - £ — — a^F^(x)-j— Z[J, J] 

m^OSi^oo il dJ(x) 6J(X) /_/_() 

1 J[dU] e~ A « detD(m q ) Z x tr[(D c + m^a^F^x)] 

= lim lim — — - — : -. — - — — (52) 

™,-o^ooO J[dU] e~ A <> detD(m q ) v ; 
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where Q = (Na) 4 is the volume of the 4D lattice, Z[J, J] is the generating 
functional in (jljl. and the Dirac, color, and flavor indices have been suppressed. 
Here the matrix valued field tensor F^ v (x) can be obtained from the four 
plaquettes surrounding x on the (/}, z>) plane, i.e., 

ga 2 F^(x) ~ ^\P^(x) + P^ u (x - jj) + P^x - z>) + P M „(x - /} - z>) 
-P^) - p t»(x - A) - PjU^ - - Plui x ~ A - #)] 

where 

P^(z) = U„(x)U v {x + /x)E/J(a: + i>)C^(z) . 



a g <q%/ q> 

0.000 i 




-0.004 - 



-0.005 1 1 1 1 1 1 

0.00 0.05 0.10 0.15 0.20 0.25 

m q a 

Figure 9: The quark-gluon condensate a^g^a^F^q) versus the bare quark 
mass m q a. The solid line is the linear fit. 

In the following, we measure the quark-gluon condensate in quenched ap- 
proximation, with one fixed site at the origin (0, 0), for 100 gauge configurations 
generated with Wilson gauge action at (3 — 6.0 on the 16 3 x 32 lattice. The 
results are plotted in Fig. |H]for 13 bare quark masses, and they are well fitted 
by the straight line 

- a 5 g(qa^F^q) = 5.59(12) x 10" 4 + 0.0150(0) x (m q a) . (53) 
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Now we take the value 5.59(12) x 10 4 at m q a = to be the quark-gluon 
condensate, then we obtain 

g(q^uF^ u q) = -0.0170(4) GeV 5 (54) 

with the inverse lattice spacing a -1 (J23j) determined in Section 2. Using the 
one loop renormalization [21], we transcribe the quark-gluon condensate (p>f)) 
to MS at scale /i = 2 GeV, 

g(qcj f , u F IIU q} MS (2 GeV) = -(434 ± 4 MeV) 5 , (55) 

Thus the ratio of the quark-gluon condensate (jB3|) to the quark condensate 
(EH) is 



MS 

(2 GeV) = 0.98(2) GeV 2 , (56) 

which becomes 0.92(2) GeV 2 at // = 0.5 GeV, in agreement with the estimate 
0.8 ± 0.2 GeV 2 [34 in the QCD sum rule. 

5 Summary and Discussions 

In this paper, we have determined several quantities [the topological suscepti- 
biltiy (jnij), the parameters 5, Ai,B (J2HJ)-(|II0J) in the pseudoscalar meson mass 
formulas, the light quark masses ()44|) - (j43j) . the quark condensate ()51|) . and 
the quark-gluon condensate (|55jl ] in quenched lattice QCD with exact chiral 
symmetry, with quark fields defined by the boundary modes 0-© of the op- 
timal DWF. The propagator of the 4D effective Dirac operator, D _1 (m g ), (fTTj) 
is solved with Neuberger's 2-pass algorithm [TT], for 100 gauge configurations 
[TH] generated by Wilson gauge action at (3 = 6.0 on the 16 3 x 32 lattice. Our 
results provide strong evidence that lattice QCD with exact chiral symmetry 
indeed realizes the quenched QCD chiral dynamics. 

The coefficient of quenched chiral logarithm (5 = 0.164 ± 0.013) extracted 
from the pion mass agrees very well with that (5 = 0.16 ± 0.02) obtained from 
the topological susceptibility. The agreement provides a consistency check of 
the theory. Also, both of them are in good agreement with the theoretical 
estimate 5 ~ 0.176 in quenched chiral perturbation theory. 

Even though the size of our lattice is insufficient for computing the light 
hadron masses at the physical (e.g., ~ 140 MeV) scale, we can determine 
the values of the parameters S, A x , B fl2^ - (|3T]|) in the pseudoscalar meson mass 
formulas of qxPT, and then use these formulas to evaluate the pseudoscalar 
masses at the physical scale, as well as to determine the light quark masses 
m Utd and m s , @3J-@5J. 



Ml 



q) 
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Using the Gell-Mann-Oakes-Renner relation, 

4m (t/S^ 



ml 



f 2 

J 7T 



and Eqs. (|37|) and ()13jl. we can estimate the chiral condensate at the zeroth 
order 



2 



CP 

(uu) = J f- = -0.00942(18) GeV 3 (57) 



m ( 2m 



(ss) n = (uu) n — — £ - 1 = -0.0108(3) GeV (58) 
u m s \ mi J 

where (J2BJ), (|HEJ), and the experimental inputs (m„ = 135 MeV, mx = 495 
MeV) have been used. Using one-loop renormalization (|43|) . we transcribe (JH7j) 
to MS at scale /i = 2 GeV, 



(mz)^ S (2 GeV) = -(223 ± 3 MeV) 3 , (59) 
(ss)^ S (2 GeV) = -(233 ± 6 MeV) 3 . (60) 

It is expected that the values of the quark condensates (JSHJ) - (jHOj) are different 
from (|51j) . since the formers are zeroth order estimates using the pseudoscalar 
mass formulas with parameters determined by the long distance behaviors 
of the pion propagator, while the latter is determined by the short distance 
bahaviors of the quark propagator. 

The quark-gluon condensate as well as the chiral condensate are the basic 
features of the QCD vacuum. In the framework of the QCD sum rules, they 
are introduced as parameters in the nonperturbative generalization of opera- 
tor product expansion, and are related to hardonic properties. Therefore, it is 
important for lattice QCD to determine these quantities from the first princi- 
ples. Our determination of the quark-gluon condensate (J55|) is the first result 
of this quantity from lattice QCD with exact chiral symmetry. The agree- 
ment of the ratio (of the quark-gluon condensate to the chiral condensate), 
M W(0.5 GeV) = 0.92(2) GeV 2 , with the estimate 0.8(2) GeV 2 in the QCD 
sum rule is encouraging, which suggests that lattice QCD with exact chiral 
symmetry does provide a viable framework to tackle the low energy physics of 
QCD pertaining to the nonperturbative QCD vacuum. 
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